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1. Introduction 

Since the discovery of quasicrystals in the early eighties |1|, |2|, 0, |j] considerable attention 
has been paid to the magnetic properties of these materials. While many quasicrystals 
contain atoms (such as Fe, Mn, or rare-earth elements) that carry local magnetic 
moments, these are usually screened very effectively, and consequently one finds a weak 
paramagnetic or diamagnetic behaviour, see e.g. H |J. Recently, however, there has 
been ample experimental evidence for magnetic ordering in quasicrystals, including 
ferrimagnetic |7j], ferromagnetic f8|, anti-ferromagnetic 0, and spin-glass behaviour 
flO| , |i~2| , though some results are still discussed controversially, see e.g. [|14], [15|], 
in particular with regard to the importance of crystalline phases present in the samples. 

Even before magnetic ordering in quasicrystals had been observed experimentally, 
theoretical investigations on the influence of quasiperiodic order on magnetic properties 
commenced. In most cases, the models considered were either one-dimensional quantum 
spin chains with aperiodic sequences of coupling constants or classical Ising models on 
two-dimensional quasiperiodic graphs; we refer the reader to the recent review [0 for 
a rather complete list of references. Recently, a symmetry classification scheme for 
magnetically ordered quasicrystals has been proposed [ |H]j . 

In this context, it is one of the central questions whether quasiperiodic order 
influences the universal properties at the phase transition, such as the critical exponents, 
in comparison to the periodic case. There is a heuristic criterion due to Luck [ITS] on 
the relevance of aperiodicity. According to this criterion, the "topological disorder" 
encountered in two-dimensional quasicrystals, generated by the cut-and-project method, 
is irrelevant; and hence an Ising model on a quasiperiodic tiling should belong to the 
same universality class as the Ising model on the square lattice. Clearly, non-universal 
properties do in general depend on the particular system under consideration. For 
instance, the location of critical points of lattice models depends in a systematic way on 
the structure of the graph on which the model is defined, see |19| and references therein. 

In this article, we consider high-temperature expansions of the free energy for zero- 



field Ising models on two planar quasiperiodic graphs, the decagonal Penrose |20|, |21| and 
the octagonal Ammann-Beenker [^, ^] tiling. The technique of high-temperature 
expansions is well known, see e.g. p5| ; it was developed several decades ago and has 



since been applied to a variety of periodic lattices in both two and three dimensions. 
With regard to previous work on high-temperature expansions of quasiperiodic Ising 
models, we are only aware of two articles by Abe and Dotera [2E[ ^71] who compute the 



expansion of the free energy up to the eighth order for the Penrose tiling and its dual, 
and of a few numerically calculated expansion coefficients for the susceptibility for the 
Penrose case Employing a systematic procedure, we are able to compute the exact 
values of the coefficients up to the 18th order for both the Penrose and the Ammann- 
Beenker tiling. This requires much more effort than the calculation for periodic lattices, 
because the number of graphs that one has to take into account grows tremendously 
with the order. Although our expansions are still not yet sufficient to extract good 
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estimates for the critical temperatures or the critical exponents, we can show that our 
results are consistent with those obtained by different methods. 

Presently, the most accurate data on the transition temperature and the critical 

Besides graphical 



31, 32 



exponents stem from Monte-Carlo simulations |29|, |3T| 
expansions and Monte-Carlo simulations, further methods have been employed to gain 
information about the critical behaviour of quasiperiodic Ising models. First of all, 
exactly solvable cases can be constructed as, for instance, the Ising model on the 
so-called labyrinth tiling [ 33| , see also jl6[ |34| for further examples. These models 



correspond to particular choices of coupling constants, restricted by the requirement of 
integrability, and thus might not be representative for the general situation. For the 
solvable models based on the idea of "Z-invariance" , see |l6j and references therein, the 
critical behaviour necessarily is the same as for the periodic case, but one does not get 
a clue whether this extends to the general case, or whether it is at least the generic 
situation. Secondly, there is an interesting approach using Lee- Yang zeros f3~5fl , which 
are complex roots of the partition function in certain variables. Simon and Baake [36[ 
calculated the zeros of the partition function for a large patch of the Ammann-Beenker 
tiling numerically and drew conclusions about the critical temperature and the critical 
exponents. Furthermore, renormalization group techniques were applied to study the 
Ising model on two-dimensional quasiperiodic tilings [P5| . In that case, one exploits the 
self-similarity of quasiperiodic tilings, which translates into a renormalization procedure 
that, however, can only be treated approximately in general. We note that for one- 
dimensional quantum Ising chains with aperiodically modulated coupling constants, 
corresponding to two-dimensional layered Ising models, renormalization techniques may 
yield exact results for the critical behaviour |39], |(J. In this case, the modulation is 
one-dimensional, and in accordance with Luck's criterion |TSf one finds that the critical 
behaviour depends on the fluctuations of the aperiodic sequence of coupling constants. 

So far, all results appear to be in accordance with Luck's criterion, including a 
recent Monte-Carlo study of the three-state Potts model on the Ammann-Beenker tiling 
Still, most approaches are based on numerical or approximative treatments 
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In 

contrast, it is our aim to obtain exact values of the coefficients for the high-temperature 
expansions in the present paper, which is organized as follows. In the subsequent section, 
we briefly recall the graphical high-temperature expansion of Ising models. In section ^, 
we discuss the generation of subgraphs of quasiperiodic tilings and the computation of 
their occurrence frequencies. Then, in section |J we present our results for the coefficients 
of the high-temperature expansions for the Penrose and the Ammann-Beenker tiling. 
The corresponding implications for the critical behaviour are discussed in section [5| In 
section |, we compare our results with exact calculations of the partition functions of 
periodic approximants. Finally, in section [7|, we present our conclusions. 
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2. High-temperature expansion 

We now give a brief account of the high-temperature expansion for the free energy of 
an Ising model on a graph without an external field |[25|| . Let us consider a finite graph 
Q containing N sites (vertices) with M neighbour pairs of vertices connected by bonds. 
We emphasize that, throughout the paper, the notion of neighbouring vertices refers to 
vertices connected by a bond, and not to the geometric distance between the vertices. 
For instance, in the Penrose tiling discussed below, the short diagonal of the small 
rhomb corresponds to the smallest distance between vertices, but does not constitute 
a bond. At a vertex j, we place an Ising spin o~j G {±1}; and two spins o~j and at 
located at neighbouring vertices j and k interact with a coupling constant J which we 
assume to be independent of the position. Hence the energy of a spin configuration 
a — {<7i, cr 2 , . . . , a N } on Q is given by 

E(<t) = -Jj2°j°k (2.1) 

where we sum over all neighbour pairs (j, k) connected by bonds as mentioned above. 
The logarithm of the partition function 

Z(G)= ^Texph/^)] 

= [cosh(/3J)] M II t 1 + °i° k tanh (/3 J)] (2.2) 
is, apart from a factor —1/(3, the free energy. It can be expanded as 

oo 

\nZ(Q) = N ln2 + M ln[cosh ((3 J)} + N 9n w n , (2-3) 

n=l 

where (3 = 1/ksT with Boltzmann's constant ks and temperature T. The expansion 
variable 

w = tanh(/3J) , (2.4) 

is small for high temperature, hence the notion high-temperature expansion. The 
expansion coefficients g n are related to the number of subgraphs of Q containing n 
bonds. 

The terms in the expansion ( |2.3j ) can be rearranged in a different fashion which is 
more convenient for our needs (see page 382 in P5| ) 

In Z{Q) = In Z{Q) - N In 2 - M In [cosh {(3J)\ 

oo 

= N^2g n w n = ^2(cr;g)kr{w), (2.5) 

n=l r 

where we now sum over all connected subgraphs c r of Q. The quantity (c r ;G) denotes 
the so-called lattice constant of c r mQ, counting the number of ways c r can be embedded 
in Q. The weight functions k r {w) depend only on c r , not on Q. In our case, without 
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external field, we can restrict the sum to so-called star graphs. These are graphs that 
cannot be dissected into two disjoint subgraphs by eliminating a single vertex. 

The weight functions k r (w) in equation (j2.5|) can be calculated from the partition 
function Z(c r ) of the subgraph c r . For this aim, let us generate all star subgraphs and 
arrange them in a sequence {c r } r= \^,... such that c s cannot be embedded in c r for r < s. 
In other words, the lattice constant (c s ; c r ) may be non-zero only if s < r, which, in 
general, does not determine the sequence uniquely Having arranged the subgraphs in 
such a way, the expansion (|2.5| ) for a subgraph c r gives 



\nZ(c r ) = ^(c s ; Cr) k s (w) (2.6) 



s=l 



and, taking into account that (c r ; c r ) = 1, we obtain the corresponding weight k r {w) 



r-l 



k T {w) = \n.Z{c r ) — ^^(c s ; c r ) k s (w) (2.7) 

s=l 

expressed in terms of lattice constants (c s ; c r ) and weights k s (w) with s < r. Therefore, 
we can compute the weights k r (w) successively provided we know the partition function 
Z(c r ) and the lattice constants (c s ; c r ) of all star graphs c s that are subgraphs of c r . 
We note that we can rearrange the sum in equation ( |2.5|) as 

oo 

^z(g) = EEE^^SM (2-8) 

n=3 r s 

where r labels closed loops lr consisting of n bonds, and are all possible complete 
"fillings" of the loop l r ■ By "fillings" of a loop we mean all proper subgraphs of Q 
which have the loop as their boundary. Here, the functions k^) (w) have the form 

kl n }(w) = w n + 0(w n+1 ) , (2.9) 



hence truncating the sum over n in equation ( [2.8|) yields all terms in the expansion up 
to nth order in w. The calculation of the weight functions k r n ](w) can be performed in 



analogy to that of the weight functions k r (w) ( |2.7|) . 

In summary, in order to calculate the high-temperature expansion (|2.8| ) of the Ising 
model to order n max we have to perform the following steps: 

generate all loops lr in the graph Q consisting of n < n mSLX bonds; 
construct all fillings Cr™2 of l r n ^; 

calculate In Z(ci™]), the logarithm of the partition function for the subgraphs c^; 



(ii 
(iii 

(iv 

(v 
(vi 



calculate the lattice constants (c r J;Q) and {c r , s ,; cf )S ;); 
compute the weight functions ki 1 ^) (w) by successive use of the analogue of ( |2.7D 
calculate the expansion (|2] 



We are now in the position to apply this scheme to the case of quasiperiodic graphs. 
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3. Frequencies of subgraphs of quasiperiodic tilings 



In fact, we want to obtain the expansion fl2.5|) for the Ising model on an infinite 
quasiperiodic graph Q. Therefore, we have to compute the corresponding "averaged 
lattice constants" per vertex 

(cr',g) := lim — (c r ;G N ) (3.1) 

N—KXl 1\ 

where Qn denotes finite patches with N vertices approaching the infinite graph Q. In 
other words, we need to calculate the occurrence frequency of a subgraph c r in the 
infinite graph Q. The main challenge now is to compute these quantities for a given 
quasiperiodic graph and all of its subgraphs up to a certain size. 

For quasiperiodic graphs generated by the cut-and-project method |42|] the 



frequencies of subgraphs can be computed exactly. In the cut-and-project method, 
one starts from a higher- dimensional periodic lattice, and projects a certain part of 
it onto a lower- dimensional "physical" or "parallel" space E\\. For the two cases of 
interest, the Penrose and the octagonal Ammann-Beenker tiling, the lattices have to be 
at least four-dimensional, the minimal choice being the root lattice A 4 for the Penrose 
case [J3| and the hypercubic lattice Z 4 for the octagonal case p3 |. The root lattice 



A 4 can be considered as a sublattice of Z 5 , wherefore the latter, albeit not minimal, 
is frequently used to generate the Penrose tiling. The physical space E\\ is determined 
as an invariant subspace with respect to the relevant subgroup (in our examples the 
dihedral groups D 5 and D$, respectively) of the point group of the periodic lattice. Its 
orthogonal complement, the perpendicular space E±, is then also an invariant subspace 
of this symmetry. The quasiperiodic tiling is now obtained by projecting all those lattice 
points onto En whose projection onto E± falls into a certain set called the "window" 
or "acceptance domain" A. In the minimal case, this acceptance domain has the same 
dimension as E±; however, if we project the Penrose tiling from the hypercubic lattice 
Z 5 , the perpendicular space is three-dimensional and the acceptance domain consists 
of four regular pentagons P m (m = 1,2,3,4) situated on equidistant, parallel planes, 
and two isolated points (Pq and P5), see figure [l| For the Ammann-Beenker tiling, the 
situation is simpler; the acceptance domain, which is obtained as the projection of the 
four-dimensional hypercube to E±, is a regular octagon O. 

Now, considering an arbitrary motive c consisting of a collection of p points 
c = {Ty ; 1 < i < p} in physical space, we can compute its occurrence frequency, 
i.e., how often translated copies of the point set occur in the infinite tiling. Associated 
to the set c of points in physical space is a corresponding acceptance domain A(c) C A 
in perpendicular space, obtained by intersecting p copies of the acceptance domain A 
shifted appropriately with respect to each other. This corresponds to the acceptance 
domain filled by choosing a reference point of the motive c, and, for all occurrences 
of the motive in an infinite tiling, lifting the positions of this reference point to the 
higher- dimensional lattice and projecting to E±. Hence, the area of A(c), divided by 
the area of A, is the occurrence frequency of our motive, as follows from the uniform 




Figure 1. The acceptance domain of the Penrose tiling, consisting of four regular 
pentagons Pi, P2, P3, P4, and two isolated points Po, P5, situated on equidistant 
parallel planes in the three-dimensional space E±. The polytope spanned by the lines 
is the projection of the five-dimensional hypercube to E±. 



distribution on the acceptance domain, see H] and references therein. 

In the Penrose case, the acceptance domain A(c) consists of four pieces A m (c) C P„ 
(m = 1, 2, 3, 4) which have to be taken into account. They are given by 



A m ({rf}) = p|{p f 



(3.2) 



where P m = if m ^ {0, 1, 2, 3, 4, 5}. The coordinates rj| G E\\ and G E± have the 
form 

/ 1 \ 



E 

i=o 



n 



(0 



cos 
sin 



27TJ 
27Tj j J 



,(*) 



E 

i=o 



(0 



cos^f 
\srn¥/ 



(3.3) 



with integer coefficients n * which correspond to the coordinates of the lattice point in 



Z 5 that projects to . The first component of r\_ , 

j=0 



J 



n 



(3.4) 



denotes the so-called translation class of the point ri , which just labels the part of 
the acceptance domain P t (i) where the corresponding perpendicular projection lies. In 
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Figure 2. The "fattest" loop of length 8 in the Penrose lattice and its acceptance 
domain (black polygons) with respect to the reference point marked by the circle (o). 
The area fraction is t — 8/5 ~ 0.0180 where t = (1 + \/5)/2 is the golden ratio. The 
symmetry factors read R = 5 and S = 1 , thus the occurrence frequency of this loop in 
the Penrose tiling, in an arbitrary orientation, is 5t — 8 ~ 0.0902. 




Figure 3. The same as figure |^, now for the "fattest" loop of length 10. Here, the 
area fraction is (14r — 22)/5 ~ 0.1305, and the symmetry factors read R = S = 1. 



figures |2] and || we show two examples where the motives are the "fattest" loops, in 
terms of the enclosed area, of length 8 and 10 in the Penrose tiling that contribute to 
the high-temperature expansion. 

For the eight-fold Ammann-Beenker case there is only one acceptance domain O, 
hence 

A({rf}) = C\{°- r ±} > ( 3 - 5 ) 

i 

where the projections to En and E± are given by 



r 



(0 _ ( cos f ^ _ STAi) ( cos — 



Here, G Z denote the coordinates of the lattice point in Z 4 that projects to ri . 

The acceptance domains of a motive c are intersections of convex polygons and 
hence themselves polygonal, see figures |2| and |3|. It is readily seen that the coordinates 
of the vertices of the acceptance domains belong to certain extensions of the field of 
rational numbers Q. For the Penrose tiling, one has to perform the calculation in the 
field 



r. 



V2 + r) = | a + h\Jl + r + cr + dr\/2 + r a, 6, c, d G Q} (3.7) 
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(-12+5A,)/8 




(22-9X)/8 




(22-9X)/S 




Figure 4. A loop of length 10 in the Ammann-Beenker tiling that can be filled in 
three different ways. The corresponding occurrence frequencies of the filled patches, 
obtained from the area fraction of the acceptance domains shown on the right, are 
given below the patches. They add up to 4 — 13A/8 ~ 0.0769 which is the frequency 
of the (empty) loop in the Ammann-Beenker tiling. The encircled node denotes the 
reference point. 




(-12+5X)/8 (22-9X)/8 (22-9X)/8 (-12+5A,)/8 (-12+5X)/8 



D 



Figure 5. Same as figure ^ for another, reflection-symmetric loop of length 10 which 
can be filled in five ways obtaining three reflection-symmetric patches and one pair of 
patches that map onto each other under reflection. Here, the frequency of the (empty) 
loop is 1 - 3A/8 ~ 0.0947. 



where r = (1 + \/h)/2 is the golden ratio, satisfying the quadratic equation r 2 = r + 1. 
For the Ammann-Beenker case, the corresponding number field is 

Q(A) = {a + b\ | a,b E Q} (3.8) 

where A = 1 + y/2 is the "silver mean" that is a solution of the quadratic equation 
A 2 = 2 A + 1. Therefore, in order to compute the occurrence frequency of a given motive 
c in the tiling Q, we have to determine the area of the acceptance domain carrying out the 
calculation in the appropriate number field. The averaged lattice constant (c ; Q) is the 
occurrence frequency of c summed over all possible orientations of the motive. In these 
quasiperiodic tilings, the frequencies of motives are independent of their orientation, 
hence we do not need to calculate them separately, but just have to count how many 
orientations of the motive occur in the tiling. 

Let us focus on the Penrose tiling as an example. Rotating the motive c by an angle 
nk/b (k e Z) essentially corresponds to a rotation of the acceptance domain by 2nk/5. 
Furthermore, also the mirror image c of the motive c occurs with the same frequency, 
since the corresponding acceptance domains A m (c) are just — y4 5 _ m (c). Therefore, in our 
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expansion ( |2.5| ), it is advantageous to jointly consider graphs which are mirror images 
of each other because they give the same contribution. For this reason, we assign two 
symmetry factors R G {1,2,5, 10} and S G {1,2} to each of these graphs, R counting 
the number of rotations by angles irk/5 which do not map the graph onto itself, and 
S — 2 if reflection does not map the graph onto itself or onto a rotated copy of itself, 
compare figures fj and ^. The averaged lattice constant (c;Q), as defined above, is thus 
R times the area fraction obtained for a fixed orientation of the graph c. Multiplying 
(c ; Q) by the factor S, we can restrict the sum in equation (|2~5| ) to graphs that are 
non-equivalent under reflection. 

Eventually, we have to consider all star subgraphs of the quasiperiodic tiling, 
corresponding to all possible fillings of loops. In contrast to the case of simple planar 
(periodic) lattices, a loop in the quasiperiodic tilings can have several fillings, which may 
occur with different frequencies. In figures f| and ||, the possible fillings, together with 
the corresponding frequencies, of two exemplary loops in the Ammann-Beenker tiling 
are shown. In orthogonal space, the different fillings correspond to a dissection of the 
acceptance domain of the loop into non-overlapping parts, see figures f| and |5]. 

In order to avoid confusion, we would like to point out once more how our 
frequencies are normalized, i.e., what the numbers given in figures (| and [| really mean. 
We emphasize that the frequency we compute is not the frequency of a particular loop 
of length n among all loops of the same length. Instead, it gives the probability that a 
randomly chosen vertex belongs to the particular loop, in an arbitrary orientation. 



4. Expansion coefficients for the Penrose and the Ammann-Beenker tiling 

The Penrose and the Ammann-Beenker tiling are both bipartite graphs, which means 
that all closed loops have an even number of edges, and at least four. Therefore, for 
zero magnetic field, only even powers of w occur in the expansion ( |2.5| ) that takes the 
form 

1 oo 

F(w) = lim -\nZ(G N ) = Y> 2n ™ 2 " (4.1) 

N— >oo iv 4 

n=2 

where Gn denotes a finite patch of the quasiperiodic graph Q containing N vertices, 
and F(w) is, apart from a factor — 1/(3, the free energy per vertex. We calculated the 
expansion coefficients g2 n up to 18th order in w for both the Penrose and the Ammann- 
Beenker tiling. The results are presented in table [1]. 

As a by-product, we obtain information on another interesting physical model, 
namely the problem of self-avoiding polygons, or closed self-avoiding walks, on the 
quasiperiodic tiling. The quantities of interest are the sums S^n of the occurrence 
frequencies of all order-2n loops which are presented in table |2|. Here, is nothing 
but the mean number per vertex of closed self-avoiding walks with 2n steps, i.e., 
random walks with 2n steps that never return to a vertex visited before, except for 
the end point which equals their starting point. For regular and recently also for "semi- 
regular" lattices, there exist data for rather large values of n in the literature ||45| ; the 
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Table 1. The expansion coefficients g2n of the free energy of the zero-field Ising model 
on the Penrose and the Ammann-Beenker tiling. The values for the square lattice are 
included for comparison. 



2n Penrose tiling Ammann-Beenker tiling Square lattice 



4 


1 




1.00 


1 




1.00 


1 


6 


9-4r 


CO 


2.53 


A 


CO 


2.41 


2 


8 


12| -4r 


co 


6.03 


47± - 17 A 


co 


6.46 


^2 


10 


25lf - 1441 r 


co 


18.28 


138- 50 A 


CO 


17.29 


12 


12 


731§-416t 


CO 


58.73 


803| - 310§ A 


CO 


53.72 


37§ 


14 


1784 - 969 t 


co 


216.13 


-1220 + 586 A 


CO 


194.73 


130 


16 


-27821§ + 17750 r 


co 


898.35 


96| + 2951 A 


CO 


810.15 


490± 


18 


-124027 + 79078§ r 


co 


3924.97 


-108706 + 46566± A 


CO 


3715.07 


1958§ 



Table 2. The mean number (per vertex) of self- avoiding 2n-step polygons S271 on the 
Penrose and the Ammann-Beenker tiling, and on the square lattice. 



2n Penrose tiling Ammann-Beenker tiling Square lattice 



4 


1 = 


1.00 




1 = 


1.00 


1 


6 


9-4r ~ 


2.53 




A ~ 


2.41 


2 


8 


15 -4t ~ 


8.53 


50 


- 17A ~ 


8.96 


7 


10 


309| - 168±r ~ 


37.45 


142 


-44 A ~ 


35.77 


28 


12 


1066- 552 t ~ 


172.85 


1173- 


-416 A ~ 


168.69 


124 


14 


6400 - 3405 t ~ 


890.59 


1704- 


- 353 A ~ 


851.78 


588 


1G 


5093- 170 t ~ 


4817.93 


27175- 


9356 A ~ 


4587.62 


2938 


18 


75115- 29655 r ~ 


27132.20 


5992 + 


8178 A ~ 


25735.44 


15268 



square lattice numbers are series M1780 in f46fl. A related problem, the enumeration 



of self- avoiding walks on quasiperiodic tilings, was already investigated by Briggs |47 
However, his results are based on counting walks emanating from a fixed starting point, 
whereas we compute the exact average over all possible starting points for the self- 
avoiding polygons. Note that the number of walks does depend on the initial vertex; 
however, the asymptotic behaviour should be independent of this choice. 

The coefficients gm and Sm listed in tables [l| and belong to degree-2 extensions of 
the field of rational numbers, namely Q(r) for the Penrose and Q(A) for the Ammann- 
Beenker tiling, respectively. We note that for the Penrose case the frequencies of 
subgraphs, and thus the coefficients g 2n and S 2n , belong to the field Q(r), whereas 
the areas of their acceptance domains in general are elements of Q(r, \/2 + r). 

The limitation of our calculations was caused by a strong, exponential growth of 
the number of graphs which have to be taken into account. For the Penrose tiling, 
we have — even after identifying graphs that are equivalent by rotation or reflection 
- to deal with more than 300 000 different graphs contributing to the 18th order, see 
table [| and their quantity grows approximately by a factor between 6 and 7 when 
increasing the order by 2. The corresponding number of graphs for the square lattice, 
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included in table [3|, are much smaller; the sequence of these numbers is apparently not 
contained in PBfl . We generated the order-2n loops as boundaries of patches that are 
constructed iteratively by successively attaching rhombi to their surface, terminating 
the process when attaching further rhombi does not lead to new order-2n loops. By this 
procedure, we make sure that all graphs are found. However, we have to pay the price 
that topologically identical graphs are obtained repeatedly and have to be rejected, thus 
slowing down the procedure substantially. 



Table 3. The number of symmetry-inequivalent closed loops of order In contributing 
to the high-temperature expansion and the number of patches obtained by filling the 
loops. 



2rc 


Penrose 


Ammann- 


Beenker 


Square lattice 


empty 


filled 


cmpty 


filled 


empty/filled 


4 


2 


2 


2 


2 


1 


G 


6 


G 


4 


4 


1 


8 


24 


28 


17 


20 


3 


10 


143 


174 


77 


112 


G 


12 


839 


1034 


479 


743 


25 


14 


5634 


6957 


3007 


4981 


86 


1G 


37677 


46712 


20175 


35063 


414 


18 


255658 


317028 


139146 


244638 


1975 



5. Critical behaviour 



In many cases, high-temperature expansions yield good estimates of the critical 
temperature and the critical exponent of the free energy. The simplest approach, which 
is commonly used for this purpose, uses the ratio of two successive coefficients g2n/92n-2 
in the expansion | |2~5| . Assuming that the free energy F{w) behaves in the vicinity of 
the critical point w c as 

F{w) ~ (l-™ 2 /<f , (5-1) 
one can easily see from (|4.1|) that 



92n 
<?2n-2 



1 



w,. 



1 - 



1 



n 



0(n 



(5.2) 



In other words, for sufficiently large values of n, the ratios g2nl Qin-2 should lie on a 
straight line when plotted as a function of n~ x . The slope of this line and its displacement 
from the origin determine the critical point w c and the exponent k. 
We may estimate the critical temperature from the sequence 



g(2n) 



n- 



92r. 



-in -I 



9ln-2 



t 9ln-2 
92n-4 



(5.3) 



that approaches w\ in the limit n — > oo. In table ||, we show the results for g(2n) for the 
two quasiperiodic tilings under consideration and compare these with the estimates of 
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the critical point from Monte-Carlo simulations |3T], |32|. The corresponding values 
for the square lattice are included for comparison. 



Table 4. Estimates of the critical point of the Ising model on the Penrose tiling and 
the Ammann-Beenker tiling, and on the square lattice. 

g{2n) 



Penrose 



Ammann-Beenker 



a After reference [29 . 
b After reference [31 . 
c After reference [32 . 

d This corresponds to the exact value w c = \/2 — 1 



Square 



In 


tiling 


tiling 


lattice 


8 


0.5116 


0.2892 


0.3333 


10 


0.1778 


0.3725 


0.2308 


12 


0.2430 


0.1902 


0.1875 


14 


0.1543 


0.1486 


0.1752 


1G 


0.1334 


0.1264 


0.1726 


18 


0.1648 


0.1252 


0.1728 


9 

wi 


0.1563(5) a 
0.1552(6) b 


0.1566(5) c 


0.1716 d 



As one can see, the convergence of g(2n) is rather poor for the quasiperiodic tilings. 
In general, the rate of convergence is determined by additional singularities w' c G C of 
F(w) lying close to w c in the complex plane. These give a correction to g 2n / '<?2n-2 which 
behaves like 0[(w' c /w c ) 2n ] [25[. The influence of these corrections must be substantial 
in our case rendering the method rather inapplicable for us. We will come back to 
this point in section |6| below when we discuss the corresponding quantities for periodic 
approximants, compare also figure [8] that contains a plot of the ratios g2n/ 92n-2 for the 
case of the Penrose tiling. 

There is, however, another method which is more suitable for us to examine the 
critical behaviour. Let us consider a sequence of partial sums F m of the expansion (|4.1| ) 
at the critical point w c 



E 



g2nW 2 c n 



(5.4) 



If the function F(w) behaves like ( p.l|) , then the asymptotic behaviour of the coefficient 



92n 



w 



2n 



g2n of its expansion in the variable w jw c is given by #2n 



~ n 



-K-l 



for n 



00 



5lj . Therefore, for large m, we have 



E 

n=m+l 
00 



n=m+l 



92nW, 



2n 



E 

n=m+l 



92n 



n 



bm 



(5.5) 
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Figure 6. The dependence of the partial sums F m of equation (5.4) on to~ 2 for 
the Penrose tiling (O ), the Ammann-Beenker tiling (0), and the square lattice (□), 
respectively. The straight lines are least-square fits to the data, disregarding the three 
points with smallest m values. 



where b is a parameter and the last relation is obtained by approximating the sum by an 
integral. Therefore, for sufficiently large m, the values F m should lie on a straight line 
when plotted versus m~ K . In figure || we plot the partial sums F m for the Penrose and 
the Ammann-Beenker tiling, taking k = 2 and w c equal to the Monte-Carlo estimates of 
3l"l , |32fl , see also table |j. For comparison, we also included corresponding data for the 
square lattice where the exact solution is known. Apparently, the data points lie close to 
a straight line for all three cases, and the fluctuations in the data for the quasiperiodic 
tilings are not visibly larger than those for the square lattice. 

Concerning the numbers of self-avoiding polygons SVn shown in table 0, one 
considers their generating function 



G(x) = J2 S '- 



2nX 



2n 



(5.6) 



which has a critical point x c that is characterized by a cusp-like singularity; i.e., in the 
vicinity of x c one has 

, 2-a 



G(x) ~ A{x) +B(x) (1 -x 2 /x 2 c )' 



(5.7) 
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15 




Figure 7. Ratios of the numbers of self-avoiding polygons (per vertex) on the Penrose 
tiling (O ) ; the Ammann-Beenker tiling (0), and the square lattice (□), respectively. 
Square lattice data for 2n < 56 are taken from [E2I. The straight lines are obtained 



from equation (5.7), in analogy to equation (5.2); using the critical exponent a = 1/2 
and the approximate values of the critical point x c given in Wn and M5J. 



with a critical exponent a, and A(x) and B(x) are non-singular at x — x c . We note that 
the only exact result for the related problem of self-avoiding walks in two dimensions is 
obtained by the Coulomb gas approach |53| and gives a critical point x 2 c = l/(2+v / 2) and 



critical exponents a = 1/2, 7 = 43/32 = 1.34375 and v = 3/4 for the hexagonal lattice 
Frequently, the so-called connective constant fi — l/x c is given instead of x c . In [47 



estimates of the critical point x c for self-avoiding walks, which coincides with the value 
for self-avoiding polygons, are given based on enumerations of walks of at most 20 and 16 
steps for the Penrose and the Ammann-Beenker tiling, respectively. The corresponding 
critical exponent in this case is 7, and all results support the conjecture that the self- 
avoiding walk problems on two-dimensional lattices and quasiperiodic tilings belong to 
the same universality class. 

In figure [7], we show the ratios of successive numbers S^/ 5 , 2n-2 as a function of 
1/n, which, by the same arguments that led to equation (|5.2|), should lie on a straight 
line for large n. Clearly, this is true for the square lattice, whereas the data for the 
Penrose and the Ammann-Beenker tiling still show sizable fluctuations. The straight 
lines in figure [?| are the functions [1 — 5/(2n)]/x^, compare ( [5.2| ), where we used the 



critical exponent a = 1/2 and the value /i = 2.618158 53 cited in [IS] for the square 
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lattice connective constant, and the estimates x c = 0.363 and x c = 0.361 H7] for the 
critical points on the Penrose and the Ammann-Beenker tiling, respectively. Given the 



rather short sequence at our disposal, and the uncertainty in the estimates [|7 , the 
agreement for the quasiperiodic cases is reasonable, thus supporting the conjecture that 
the critical point of self-avoiding polygons on such quasiperiodic tilings is described by 
the same critical exponents as for the hexagonal lattice 

6. Partition functions of periodic approximants 

One may pose the question whether one can calculate the expansion coefficients g2 n in 
equation (|4.1| ) by a different method, thus verifying our results. Perhaps it might even 
be possible to calculate the partition function Z(Q) on certain quasiperiodic tilings Q 
exactly. Although this may seem hopeless, there exist methods to tackle this problem, 
which at least allow us to compute the partition function of general periodic lattices 
explicitly, thus also for periodic approximants of the quasiperiodic tilings. Let us briefly 
recall some exact results on partition functions on two-dimensional lattices. 

The first solution of the two-dimensional zero-field Ising model for the square lattice 
had been found by Onsager and Kaufman |05|, |54]] in 1944. Several years later, Kac and 
Ward |55[ developed a combinatorial approach in which the problem was reduced to 
the calculation of a determinant of a certain matrix K (see below) which depends on 
the lattice and the coupling constants between the spins. Although this approach was 
not rigorous, it appeared extremely plausible and it initiated numerous attempts to 
generalize this result to other lattices |5(| ^7|, [58|, [59| . Recently, Dolbilin et al [|(| 
proved the long-known formula 

Z{gf = det(i^) (6.1) 

for a zero-field Ising model on an arbitrary planar graph Q with arbitrary (in general 
site-dependent) spin coupling constants. The matrix elements K(ei, ej) of the 2M x 2M 
matrix K are labeled by oriented edges ej and e,,-, 1 < i, j < 2M. They are defined as 



K(ei,ej] 



1 if e, = ej 



-Wiexp [^(e^e})] if /(ej) = b{ej) and /(e,-) 7^ 6(e<) (6.2) 
otherwise 



where u>j = tanh(/3Jj) and Jj is the spin coupling constant along the edge ej, which is 
independent of the edge in our case, Jj = J. Furthermore, (e», a,-) denotes the angle 
between edges e^ and ej, and &(ej) and /(ej) are the starting point and the end point 
of the edge e,, respectively. If Q is periodic, the matrix K is cyclic and the determinant 
can be calculated exactly in the thermodynamic limit M — > 00. We can therefore 
apply ( |6.2[ ) and calculate Z(Q) exactly for periodic approximants of the Penrose and 
the Ammann-Beenker tiling. 

Let us now briefly describe how to generate periodic approximants of quasiperiodic 
tilings in the framework of the cut-and-project method discussed in section |[ The 
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Table 5. Expansion coefficients of the free energy for the Ising model on the Penrose 
tiling and its periodic approximants m = 1, 2, 3, 4, 5, 6 with Af vertices in a unit cell. 



m 


N 


2M 


94 


96 


98 


1 


11 


44 


1 


25/11 ~ 2.2727 


127/22 ~ 5.7727 


2 


29 


116 


1 


73/29 ~ 2.5172 


349/58 ~ 6.0172 


3 


76 


304 


1 


5/2 ~ 2.5000 


227/38 ~ 5.9737 


4 


199 


796 


1 


503/199 ~ 2.5276 


2399/398 ~ 6.0276 


5 


521 


2084 


1 


1315/521 ~ 2.5240 


6273/1042 ~ 6.0202 


6 


1364 


5456 


1 


862/341 ~ 2.5279 


4111/682 ~ 6.0279 


oo 


oo 


oo 


1 


9 - 4r ~ 2.5279 


12± - 4r ~ 6.0279 




m 


9io 






912 


914 


1 


175/11 


~ 15.909 




3145/66 ~ 47.652 


1812/11 ~ 164.73 


2 


504/29 


~ 17.379 




341/6 ~ 56.833 


6011/29 ~ 207.27 


3 


679/38 


~ 17.868 




6629/114 ~ 58.149 


16123/76 ~ 212.14 


4 


3624/199 


~ 18.211 




69833/1194 ~ 58.487 


42552/199 ~ 213.83 


5 


9496/521 


~ 18.226 




610783/10420 ~ 58.616 


112451/521 ~ 215.84 


6 


24921/1364 


~ 18.271 




160129/2728 ~ 58.698 


294347/1364 ~ 215.80 


oo 


25lf - 144±r 


~ 18.279 




73lf -416r ~ 58.731 


1784- 969r ~ 216.13 



m 


916 






918 






1 


29439/44 




669.07 


95119/33 




2882.4 


2 


100769/116 


rv_i 


868.70 


342484/87 


rv_i 


3936.6 


3 


33325/38 


rsj 


876.97 


222817/57 


rsj 


3909.1 


4 


709087/796 


rsj 


890.81 


2345981/597 


r^j 


3929.6 


5 


1867989/2084 


rsj 


896.35 


6128605/1563 


r\j 


3921.1 


6 


1223683/1364 


r^i 


897.13 


4015369/1023 


r^j 


3925.1 


oo 


-27821§ + 17750r 




898.35 


-124027 + 79078§r 




3925.0 



acceptance domain A and the projection onto perpendicular space E± are altered 
in a way that corresponds to replacing the irrational numbers r and A by rational 
approximants r m and X m . Here, for the Penrose tiling we use r m = f m +i/ f m where 
f m+1 = f m + fm-!, and / = 0, /i = 1 are the Fibonacci numbers, and lim m ^ 00 r m = r. 
Analogously, one defines rational approximants X m = g m+ i/g m with the "octonacci 
numbers" g m+ \ = 2g m + g m -\ and (? — 0, g\ — 1, and Xvaim^^ \ m = A for the case of 
the Ammann-Beenker tiling. 

In this way, one obtains periodic approximants of the Penrose tiling with unit cells 
containing M = 11, 29, 76, 199, 521, 1364 vertices for m— 1,2, 3, 4, 5, 6, respectively. The 
unit cells of the periodic approximants of the Ammann-Beenker tiling with m— 1, 2, 3, 4 
contain Af = 7,41,239,1393 vertices. For both tilings, the number of oriented edges 
is 2M = 4jV, because each vertex has exactly four neighbours. We note that the 
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Table 6. Same as table [|, but for approximants m = 1, 2, 3, 4 of the Ammann-Bcenkcr 
tiling. 



m 


Af 




2M 


54 


56 


58 


1 


7 




28 


1 


17/7 ~ 2.4286 


87/14 ~ 6.2143 


2 


41 




164 


1 


99/41 ~ 2.4146 


529/82 ~ 6.4512 


3 


239 




956 


1 


577/239 ~ 2.4142 


3087/478 ~ 6.4582 


4 


1393 




5572 


1 


3363/1393 ~ 2.4142 


17993/2786 ~ 6.4584 


oo 


oo 




oo 


1 


A ~ 2.4142 


47± - 17A ~ 6.4584 




m 




510 








514 


1 




115/7 


~ 16.429 




2189/42 ~ 52.119 


1395/7 ~ 199.29 


2 


708/41 


~ 17.268 




13183/246 ~ 53.589 


7994/41 ~ 194.97 


3 


4132/239 


~ 17.289 




77029/1434 ~ 53.716 


46542/239 ~ 194.74 


4 


24084/1393 


~ 17.289 




74832/1393 ~ 53.720 


271258/1393 ~ 194.73 


oo 


138 


- 50A 


~ 17.289 




803± - 310±A ~ 53.720 


-1220 + 586A ~ 194.73 



m 


516 




518 




1 


22815/28 


~ 814.82 


78479/21 


~ 3737.1 


2 


132885/164 


~ 810.27 


153121/41 


~ 3734.7 


3 


774507/956 


~ 810.15 


2664121/717 


~ 3715.6 


4 


4514157/5572 


~ 810.15 


739303/199 


~ 3715.1 


oo 


96§ + 295±A 


~ 810.15 


-108706 + 46566^ A 


~ 3715.1 



approximant m + 1 contains about r 2 = r + 1 ~ 2.618 and A 2 = 2A + 1 ~ 5.828 as many 
vertices and bonds as the approximant m for the Penrose and the Ammann-Beenker 
case, respectively. 

For each periodic approximant, we define a matrix K labeled by oriented edges e[ 
and e'j with starting point in the unit cell. It is related to the Kac-Ward matrix K (|6.2j ) 

as 

K(e' i ,e' j ) := Kfa, e,) exp (— ifciAi) exp (-i/c 2 A 2 ) . (6.3) 

Here, we can assume that e- = ej starts in the unit cell, and e'j equals ej modulo the 
unit cell, i.e., if b(ej) = £jV + rjjW, where V and W are the base vectors spanning 
the unit cell, then b(e'j) = frac(^j)V + frac(r/j)W / where frac(x) denotes the fractional 
part of x. The integers A x and A 2 are the integer parts and [%J, respectively, 
and K depends on the "wave vectors" k\ and fc 2 which are real numbers. Due to the 
fact that the Kac-Ward matrix K of the periodic approximant is cyclic, its determinant 
can be expressed as a product of the determinant of K over all values of k\ and fc 2 , 
corresponding to a reduction to the unit cell. Calculating the logarithm of Z(Q), and 




Figure 8. The ratios <72n/<72n-2 of expansion coefficients for the square lattice (□), 
the first (+) and the second (x) periodic approximant of the Penrose tiling, and for 
the Penrose tiling (• ), respectively. Lines are meant as guides to the eye only. 



taking relation ( |6.1D into account, one obtains 

In Z(G) = ^J J lndetK{k h k 2 ) ^hdk 2 . (6.4) 
o o 

Let us now expand this equation in a series with respect to w and compare it with 
the high-temperature expansion (f4~T|). For this purpose, we exploit the fact ( |6.2| ) that 
the (finite-dimensional) matrix K(k 1 ,k 2 ) has a form K(ki,k 2 ) = 1 + wL(k 1} k 2 ) where 
L{k\i k 2 ) has zero trace. Therefore, using 

det [1 + wL{ki, k 2 )} = det exp In [1 + wL(ki, k 2 )] 

= exp tr In [1 + wL(k 1 , k 2 )] (6.5) 

one obtains, expanding the logarithm in powers of wL(ki, k 2 ), 

In det K( k x , k 2 ) = tr In [1 + wL(h , k 2 )\ 



High-temperature expansion for quasiperiodic Ising models 



20 



= Y ^ tv[P(k h k 2 )]w p , (6.6) 

U p 

where, again, only even values of p yield non-vanishing contributions to the sum. 
Comparing this result with equation (|6.4|) , we derive an expression for the coefficient 
g 2n in the expansion (|4~T|) 



92n = ~Y^J J te[L In (k 1 ,k 2 )]dk 1 dk 2 (6.7) 


for the periodic approximants. We have calculated the coefficients from ( |6.7| ) for the 
leading orders in w for both the Penrose and the Ammann-Beenker tiling. The limitation 
of the calculation was due to a rapidly growing dimension of the complex matrix 
L(ki, k 2 ), which was equal to 5456 and 5572 for our largest approximants of the Penrose 
and Ammann-Beenker tiling, respectively. The results are presented in tables [5] and 
Clearly, with increasing size of the approximant, the coefficients approach those of the 
quasiperiodic system, and the coefficients of the largest approximant are already quite 
close to those of the quasiperiodic case. 

We now consider the ratios (|5.2|) for the periodic approximants of the Penrose 
tiling. The result is shown in figure ||. Although we included terms up to order 2n = 56, 
the data for the two periodic approximants do not lie on straight lines, in contrast to 
those of the square lattice. Instead, they show large fluctuations, and apparently the 
fluctuations for the smallest approximant with 11 vertices in the unit cell turn out to 
be much larger than those for the larger approximant which contains 29 vertices. It 
would be interesting to have a better understanding of this phenomenon, perhaps an 
investigation of the complex-temperature phase diagram of the periodic approximants 
can give an explanation of this observation. Again, figure |8] also shows that the data 
for the second approximant are already rather close to that of the quasiperiodic tiling, 
and one might conclude from figure |8] that the fluctuations in the ratios of expansion 
coefficients become less with increasing size of the approximant. 



7. Conclusions 



We considered the Ising model on two planar quasiperiodic graphs, the Penrose and 
the Ammann-Beenker tiling. We calculated the leading terms of the high-temperature 
expansion of the free energy exactly, using the embedding of the quasiperiodic tilings into 
higher- dimensional periodic lattices to compute the occurrence frequencies of patterns 
in the tiling. These frequencies are expressed in terms of characteristic quadratic 
irrationalities related to ten- and eightfold rotational symmetry, the golden mean 
r = (1 + y / 5)/2 and the silver mean A = 1 + \[2 for the Penrose and the Ammann-Beenker 
tiling, respectively. 

The number of graphs that contribute to a given order in the expansion grows much 
faster with the order than for a simple periodic lattice, therefore we did not go beyond 
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the 18th order in the expansion variable w = tanh (ft J) in this work. From our expansion 
alone, it is difficult to extract information about the critical behaviour. However, using 
estimates of the critical temperature obtained by other methods to analyze our data, 
we find that our expansions are in accordance with the the conjecture that Ising models 
on such planar quasiperiodic graphs belong to the Onsager universality class. 

In order to compute the expansion coefficients, we had to construct all polygons 
on the quasiperiodic graphs with up to In = 18 edges. Thus, we obtain the average 
number of such self-avoiding polygons as a by-product of our calculation. Comparison 
with earlier results on self-avoiding walks [f|7j, based on enumerating walks that start 
from a chosen vertex in the tiling, indicates that the self-avoiding polygons on the 
examples we considered belong to the same universality class. In particular, this means 
that the corresponding critical point is described by the same exponents as for the 



hexagonal lattice which are known analytically |53 



Finally, we considered periodic approximants of the quasiperiodic tilings. For these 
it is, in principle, possible to compute the free energy of the infinite periodic system 
analytically. Here, we are only interested in the leading terms of the free energy, which 
we compare with those of the infinite quasiperiodic tiling. We find that, at least for 
the leading orders in w, the rational coefficients of the approximants converge rapidly 
towards the irrational coefficients obtained for the quasiperiodic tiling. 

In conclusion, it is doubtful whether the computational effort necessary to extend 
the expansions to higher order will result in a considerable improvement of the estimates 
of the critical properties. Instead, it might be more rewarding to consider periodic 
approximants and use methods as those outlined in section |6| and |^TJ to compute 
physical quantities, such as for instance the magnetization or correlation functions. 
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